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In this paper, we examine some basic properties of the muhiple-Sine-Gordon (MSG) 

systems, which constitute a generahzation of the celebrated sine-Gordon (SG) system. We 

start by showing how MSG systems can be viewed as a general class of periodic functions. 

I Next, periodic and step-like solutions of these systems are discussed in some details. 

Ch ■ 

^ . In particular, we study the static properties of such systems by considering slope and 

phase diagrams. We also use concepts like energy density and pressure to characterize and 
CNJ ' distinguish such solutions. We interpret these solutions as an interacting many body system, 

in which kinks and antikinks behave as extended particles. Finally, we provide a linear 
|. stability analysis of periodic solutions which indicates short wavelength solutions to be stable. 
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X. INTRODUCTION 

> 

in 

I The double-Sine-Gordon (DSG) equation which is a generahzation of the ordinary Sine-Gordon 



, variety of systems in condensed matter, quantum optics, and_particle physics 

applications include the spin dynamics of superfluid ^He 



(SG) equation has been the focus of much recent investigations [ll-UJ. It has been shown to model a 

md oarticle physics 21 • Condensed matter 

incommensurate phase transitions ^ , surface structural reconstructions 0, and domain walls flfl 

X ■ n 

$H ' and fluxon dynamics in Josephson junction [10|]. 



In quantum field theory and 
and self-induced transparency 



quantum optics, DSG applications include quark confinement [11| 



I2| ]. The internal dynamics of multiple and single DSG soliton 
configurations using molecular dynamics have been studied in [^. There have also been studies 
about kink anti-kink col 
mechanics applications 



ision processes for DSG equation [13|]. One can also point to the statistical 
1^ . and perturbation theory for this system jdst ]. Following our pervious 
study on the periodic and step-like solutions of DSG equationyj, we focus on a generalization of 
this system with the self-interaction potential (see Fig{T]) 
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FIG. 1: MSG Potential for A'' = 6. The dashed curve is for e = 10, the dotted curve is for e = 1 and the 
sohd curve is for e = 0. 



V{(j)) = l + e-cos(j)-e cos{N(p), 



(1) 



where e is a non- negative constant. This potential reduces to the ordinary SG potential in the 
limit e — )• 0. A'' is an integer with = 1 and N = 2 corresponding to the SG and DSG systems, 
respectively . In this paper, we propose to study this system along the line of our previous 

work on DSG [1] for N > 2. In Section [TI] we present some general properties of the system 
including the single kink solutions, in Section Hill we examine the = 3 case, in Section ITVl N > 4 
cases are discussed, and in Section |V] the stability of periodic solutions is examined. Section IVTl is 
devoted to a summary and conclusions. 



II. GENERAL PROPERTIES OF THE MSG SYSTEM 



The Lagrangian density of the Multiple-Sine-Gordon (MSG) system is the following: 
Cmsg = ^d''4>d^(f) - [1 + e- cos<j3 - e cos{N(j))] . 



(2) 



From this Lagrangian density, the MSG equation follows 

Ocj) = -sincj) - Nes'm{N(j)). (3) 

Note that (j) is a real scalar field, and we are using a 1-1-1 dimensional spacetime (t,x) with the 
signature (-|-,— ). The potential for this system has minima at 

Cprnin = 2mT n = 0,1,2,... (4) 
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FIG. 2: Single soliton solutions and energy densities (inset) of the MSG system for {a,)N = 3, e = 10 and 
(b)7V = 4, e = 10. 

The energy-momentum tensor of the MSG equation can be obtained by the standard relation 
which follows from the Noether's theorem and is a consequence of the translational invariance of 
the Lagrangian density: 

T^- = d^cpd'-cP - g^^'LusG. (5) 

in which g^^ = diag{l, —1) is the metric of the Minkowski spacetime. 

Like the ordinary SG system, a topological current can also be defined for the MSG system 
according to 

= ^e^'^-<^' (6) 

where e^'^ is the totally antisymmetric tensor with e''^ = 1. This current is identically conserved 
{d^J^ = 0) and the total charge of any localized, finite-energy solution is both conserved and 
quantized. 

MSG equation possesses kink and anti-kink solutions which correspond to transitions between 
the spatial boundary conditions (/)(±oo) = 2n7r. The first integral of equation ([3]) for static solutions 
reads 
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in which we have used the boundary conditions (j){zizoo) = 2mT. We therefore have 



This integral cannot be carried out analytically for general N and e. A few examples of these 
solutions for = 3 and = 4 which are obtained numerically are shown in Figlll together with 
their energy density, which is equal to the Tq component of Ty . 

Multiple Sine-Gordon systems admit soliton-like kink solutions with interesting properties. 
Here, we show that an arbitrary, periodic potential with vanishing minima at = and (j) = Itt 
can be expanded in terms of MSG potentials. This property adds to our motivation for studying 
this class of nonlinear equations. 

Let us define the MSG basis functions 



N\ 



1 + e - COS </) - e cos{N(j)). (9) 



An arbitrary, periodic function /(</>) which satisfies the boundary conditions /(O) = f{2iT) = can 
be expanded in terms of this basis: 

oo oo oo oo 

/((/.) = J2 aNF'N{(t>) = + aN-cos(t>{Y^ a^v) - e ^ aNCOs{N(t>), (10) 

N=0 N=0 N=0 N=0 

in which a^v's are the expansion coefficients. At the same time, f{(p) can also be represented by a 
cosine Fourier series: 

oo 

/(0) = J]6^cos(A^0). (11) 

N=0 

of course, /(O) = /(27r) = demands that X]?^=o ~ ^■ 

Equating the coefficients of cos(A^i?!>) in ([TO]) and pT]) leads to 

oo 

bo = (1 + e) ^ OAT - eao, 

N=0 

oo 

h = - ^ QAT - eai, 

N=0 

= -euN {N > 1). (12) 
Note that Y1n=o'^N' can be easily calculated via 



^^0 2vr(l + e) 
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As an example, consider the triangular periodic function 



(j) ^ (/) ^ vr 
2tt -6 vr ^ (A ^ 27r 



with f{(j) + 2n7r) = f{(j)). The Fourier cosine series representing this function reads 

8 8 1 

/((/>) = vr - - cos((^) - — cos(3<^) - — cos(4(/.) - .... (14) 
vr yvr 47r 

The coefficients can be easily calculated from Equations (fT2]l and (fT3]) : 



TT 



2(1 + 



vr 

1( 

ai 



16(l + e)-7r2 



27re(l + e) ' 

02 = 0, 

8 

= 9^' 

...etc (15) 

The reader may get worried about the ^ dependence as e — y oo, but this is not a serious 
problem, since this is a common factor for all oat's and may be factored out. 

III. MSG FOR N = i 

In this and the next two sections, we investigate some general properties of MSG for = 3, 
and subsequently for = 4 and A^ > 4, along the line of our previous work on DSg[i|. We focus 
on the static, time-independent solution of MSG equation. To solve this equation, we use a fourth 
order Runge-Kutta method 1^ to numerically integrate: 

^ = ^ = sincP + Nesinm)- (16) 

We use (j) = TT with various values of ^ at x = 0, as our initial conditions throughout. We study 
this system for various values of A^ and e. The solutions are characterized by the first integral of 
the static MSG equation, i.e. 

p = -Tl = \i^f-n'i>), (IT) 
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FIG. 3: (a) The slope diagram for the perfodic chain of MSG sohtons for = 3 and e = 10, with 
P = -0.2859 for the sohd curve and with P -0.8750 dotted curve and P = -2.0604 for the inset. The 
three solutions have the same average density (p) but different pressure (P). The main figure shows three 
subkink solutions while the inset shows a no-subkink solution, (b) The phase diagram(|^ vs. 0) for the 
periodic chain of MSG solitons, for = 3 and e = 10. The single loop corresponds to P = —17.5, (kink 
anti-kink), while the three lobe one corresponds to P = —0.875, (three subkinks). 



with 

dP 



Note that P is different from energy density 



P = ro° = ^(f)^ + n<^); (19) 

which changes with position x. The reason for using the letter P for this quantity is that it conforms 
with "pressure" for a one dimensional perfect fluid: 

T!^ = diag{p,-P). (20) 

Here, pressure is to be interpreted as tension, rather than force per unit area The value of P 
distinguishes two different type of solutions: (i) step-like solutions for which P > 0, which are the 
"running" solution of the mechanical analog, (ii) periodic solutions for which —2(1 + e) < P < 0, 
which are the "oscillatory" solutions of the mechanical analog {isl- Here, since we have chosen 
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FIG. 4: Energy per soliton diagram for periodic chain of MSG solitons for iV = 3 and e = 10. The low 
energy portion contains no subkinks while the high energy solitons is muhi- valued and contains three-subkink 
sohtons. 

(f) = TT for all our solutions as an initial condition, the minimum value of P for periodic solutions is 
—2 for N even, where this lower bound is e-dependent for odd N. 

We first consider the simpler, step-like solutions. These solutions simply appear as in the case of 
DSG equations. The general MSG equation exhibits step-like solutions for P > which are always 
characterized by N subkinks, as in Figl21 provided that e is in the appropriate range. However, 
as in the case of DSG equation, the periodic solutions exhibit a rich set of behavior depending 
on the value of P and e, as well as A^. A simple way to imagine these solutions is to consider 
slope diagrams vs. x) or phase diagrams vs. ip). Such a set is shown in Figl3]for = 3 
and e = 10. We note that for = 3, a kink anti-kink solution (with no subkinks) as well as a 
three-subkink solution is possible. Part(a) shows the kink anti-kink solution as an inset and two 
different three-subkinks in the main panel. Part (b) shows the phase diagram where each (subkink) 
solution is represented by a "loop" which is a helpful way of distinguishing these solutions. 

A different ( "thermodynamic" ) method to characterize solutions of such systems is to consider 
the energy diagram {E vs. L) as well as equation of state diagram (P vs. In the first case 

energy per soliton (E) is plotted as a function of inter-soliton distance (L) [l9|] , while in the second 
case, P is plotted as a function of average energy density P = j Jq pdx with / is a sufficiently large 
length (includes several solitons). These figures are shown for = 3 and e = 10 in FigU] and 
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FIG. 5: Equation of state diagram for periodic chain of MSG solitons for = 3 and e = 10. The upper 
curve corresponds to the three-subkink solutions and lower one to the no-subkink solutions. 




FIG. 6: The equation of state diagram for N=3, e — 10, for step-like solutions. The P limit corresponds 
to the single kink solution. Unlike the periodic solutions, the general shape of this curve is independent of 
N. 

FigiSl respectively. As in the case with DSG equation, the energy diagram for iV = 3 shows a 
region where E is double-valued. That is, Figjl] shows a low energy solution with no subkinks, a 
region where no solutions are observed for our initial conditions, and a high energy region which 
is double valued. It is interesting to note that in the lower branch of high energy solution, as L 
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FIG. 7: (a) The slope diagram for the periodic chain of MSG sohtons for = 4 and e = 10, with 
P = -0.0398 for the sohd curve and with P ^ -0.70395 dotted curve, with inset P = -1.18080 for the soUd 
curve and P = —1.8078 dotted curve. The four solutions have the same average density (p) but different 
pressure (P). The main figure shows four subkink solutions while the inset shows two subkink solutions, 
(b) The phase vs. 0) diagram for the periodic chain of MSG solitons for = 4 and e = 10. The two 
lobe loop corresponds to P = —1.68, while the four lobe one corresponds to P = —0.38. 



increases, kink and anti-kink separate out, while in the upper branch it is the subkinks (within a 
given kink)which separate out when L increases. We note that in this three-subkink solution the 
middle subkink never changes and it is the other two (first and last) which actually extend out in 
the upper branch. 

Figl5] shows the equation of state for = 3 and e = 10. It consists of two separate branches as 
well. The upper smaller branch is the three subkink solution and the lower longer branch is the no 
subkink solution. The upper branch is similar to that of DSG solution and the lower one is similar 
to the regular SG equation. The tension (or negative pressure) decreases in this Id chain with 
increasing density due to the attractive force between the solitons in the chain. The only region 
where tension increases with increasing density is in the lower half of the upper branch which 
indicates a different type of inter-soliton force in the chain. It is worth noting that the equation of 
state diagram for the step-like solutions for > 2 is in general no different than the N < 2 case 
previously studied [l|. Here we plot the = 3 case in Figl6j 
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FIG. 8: Energy per soliton diagram for periodic chain of MSG solitons for iV = 4 and e = 10. The upper 
curve corresponds to the four subkink solutions and the lower one to the two subkink solutions. 
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FIG. 9: The equation of state diagram for = 4 and e = 10 system. The upper curve corresponds to the 
four subkink solutions and the lower one to the two subkink solutions. 



IV. > 4 CASES 



As N increases, the MSG system gets more complicated and detailed. The main reason for this 
is the development of more subkinks. The emergence of subkinks in the periodic regime depends 
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FIG. 10: The equation of state diagram for iV = 6 and (a) e = 0.03, (b) e = 0.05, (c) e = 0.07, (d) e = 0.1, 
(e) e = 1 and (f) e = 10 system. 



on the value of e as well as P. The maximum number of possible subkinks is N. Due to the 
symmetry V{(f)) = y(— 0), the number of subkinks (n) is also "doubled", i.e. 



n 



2 4 ... up to for even 

3 5 ... up to N for N odd 



One can easily draw similar diagrams for general N as the ones draw for = 3 in the previous 
section. Here, we simply show these in FigsdlHl foi' N = 4:. Here, two and four subkink solutions 
are shown, with two subkink solutions being the low energy solutions and thus lower pressure. 

The change in the equation of state diagram as e increases shows an interesting behavior. This is 
shown in FigllOlfor = 6. It is interesting to note how the multi-valuedness arises as e is increased. 
As e increases a larger range of p values allows multiple solutions. Another interesting behavior is 
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FIG. 11: The phase vs. diagram for the periodic chain of MSG sohtons, for N = 6. The nearby 
curves correspond to e = 0, 0.03, 0.06 and 0.09 from the inner to the outer ones. The inner set of loops 
corresponds to P = —1.995, as indicated, all the way to P = —0.01995 for the other set of loops. 

the emergence of cusp points with a discontinuity in ^ (the left most section in Fig(TO] (e) and (f)) 
which is due to the existence of hyperbolic unstable point in MSG equation. This discontinuity 
dose not exist for the simple SG system. Note that these cusps separate regions of solutions with 
different number of subkinks. For example, in FigllOlff) two cusps separate the 6, 4 and 2 subkink 
solutions. Another interesting point in these solutions is the maximum density solutions where 
compressibility (x = ^) becomes zero. At these point the system becomes incompressible due 
to strong interactions of kink-antikink solutions (one can see that this corresponds to a minimum 

We are also interested to see how different subkink solutions arise as e and P are changed. We 
therefore plot the phase diagram for various values of P and look for the emergence of subkinks 
as e is increased. Such a plot for = 6 is shown in Figjlll Here, each set of loops corresponds 
to a given P value. It is seen that with a small change in e subkink solutions arise. However, the 
number of subkinks is more sensitive to change in P than e. We have observed that for a given 
value of P, increasing e from zero can sometimes change the number of subkinks in a short interval. 
This interval is typically of order of e ~ 0.1. After this brief transition, the number of subkinks 
remains constant. On the other hand the number of subkink changes with increasing P (for a given 
e) until it reaches its maximum (A^) for P — > 0. 
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FIG. 12: Fixed points of the MSG equation for (a)iV = 3 and (b)A^ = 4. 



V. STABILITY AROUND FIXED POINTS 



The fixed points of the MSG equation are obtained by equating the RHS of the equation ([3]) to 
zero: 



sin + iVe sin(iV0) = 0. 



(21) 



Let us call these fixed points which lie in the range [0, 27r] as ^n^n = 1, 2, ...). The number of fixed 
points depends on both values of and e. Figll2l shows the value and number of fixed points as e 
varies for two cases of = 3 (FigfT2k) and = 4 (FigfT2b). 

In order to investigate the stability of periodic solutions around fixed points, we use the con- 
ventional linear perturbation analysis, by inserting the expression 



(x, t) = (/>n + ^cos(A;2;) cos{wt), 



(22) 



in which A{<^ (pn) is a constant and k and uj are wave number and angular frequency, respectively, 
in to the field equation ([3|). Expanding the nonlinear sine terms and keeping only terms linear in 
A, we obtain 



cj^ = /c^ + cos (j)n + A^^e cos{N(j)„ 



(23) 



In this dispersion relation, the stability is judged by the sign of u)'^ for any chosen value of k. Let 
us write Eq.(|23p in the form 

co^ = P + fniN,e). (24) 
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FIG. 13: (a) Initial periodic perturbation around the fixed point of DSG (0 = tt) for e = 0.1 and its 
FFT(inset), (b) The evolved perturbation after 200 time steps and the corresponding FFT. 

It is obvious that if /„ > 0, uj"^ is always greater than zero and the perturbation is stable. If, on 



the other hand, /„ < 0, uj"^ is positive and we have stability, only for > —fn or A < 



2n 



In 



such a case, we expect small scale perturbations to be stable and long ones to be unstable. The 
fore- mentioned linear analysis was checked numerically, by using a finite difference algorithm 17 ] . 



We have performed numerical experiments with periodic perturbations around a fixed point as 
initial condition. In this way, we could confirm stability about the minima-energy fixed points like 
<j)n = 0. Perturbations around local maxima show a more interesting behavior. It was seen that 
short wavelength modes remain stable, while long wavelength ones start to be excited and increase 
in amplitude. FigJ13l shows the fate of a periodic perturbation around = vr in the DSG system 
(A^ = 2). The initial perturbation is shown in FigJ13b together with its FFT. Fig|13b shows the 
perturbation after integrating the dynamical equation for 200 time steps. As the FFT shows, the 
initial, small wavelength perturbation persists, while the unstable, long wavelength modes grow, 
confirming our analytical result. 
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VI. CONCLUDING REMARKS 

We have examined the periodic and step-hke solutions of the MSG system. In particular, we 
interpreted the step-like and periodic solutions as a many body system of kinks and anti-kinks 
with a nonlinear interaction between the solitons. This interpretation is established by comparing 
the asymptotic P ^ case with the single-soliton solution. In this limit, each site in the chain 
conforms with the single-soliton solution, and the energy per soliton approaches the single-soliton 
energy. As the solitons get closer, the energy per soliton becomes different and the difference is 
attributed to the interaction energy. 

We observed that as A'" increases, more and more subkinks (up to N) can develop and this leads 
to a more complicated equation of state. This diagram becomes multiple- valued, depending on 
the number of subkinks, which in turn depends on the pressure P. We examined the stability of 
periodic solutions around the fixed pomt{(f> = tt and ^ = 0), obtaining the interesting result that 
within the linear approximation, short wavelength periodic solutions oscillate like standing waves 
(signalling stability), while long wavelength periodic solutions are unstable. 



M. Peyravi, A. Montakhab, N. Riazi and A. Gharaati, Eur. Phys. J. B 72, 269-277 (2009). 

S. Burdick, M.El-Batanouny and C. R. Willis, Phys. Rev. B 34, 6575 (1986). 

K. Maki and P. Kumer, Phys. Rev. B 14, 118 (1976); 14. 3290 (1976). 

Y. Shiefman and P. Kumer, Phys. Scr. 20, 435 (1979). 

K. M. Leung, Phys. Rev. B 27, 2877 (1983). 

O. Hudak, J. Phys. Chem. 16, 2641 (1983); 16, 2659 (1983). 

M. El-Batanouny, S. Burdick, K. M. Martini and P. Stancioff, Phys. Rev. Lett. 58, 2762 (1987). 
E. Magyari, Phys. Rev. B 29, 7082 (1984). 

J. Pouget and G. A. Maugin, Phys. Rev. B 30, 5306 (1984); 31, 4633(1984). 

N. Hatakenaka, H. Takayanagi, Y.Kasai and S. Tanda, Physica B 284-288 (2000) 563-564. 

T. Uchiyama, Phys. Rev. D 14, 3520 (1976). 

S. Duckworth, R. K. Bullough, P. J. Caudrey and J. D. Gibbon, Phys. Lett. 57 A, 19 (1976). 

V. A. Gani and A. E. Kudryavtsev, Phys. Rev. E 60, 3305 - 3309 (1999). 

M. Croitoru, J. Phys. A: Math. Gen. 22, 845-863 (1989). 

C. A. Popov, Wave Motion. 42(1), 309-350 (2006). 

N. Riazi and A. R. Gharaati, Int. J. Theor. Phys. 37, 1081 (1998). 

E. Kreyszig, Advanced Engineering Mathematics, John Wiley and Sons, NewYork(1983). 

This problem has a simple mechanical analog if the following change of variables are considred x - 



16 



(j) — > X, for a particle of unit mass. 
[19] Here, we use the word soliton interchangably with kink. The distance between successive kinks (L) is 
well-defined, and therefore energy per soliton {E) is well-defined as well. 



